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ABSTRACT 


The primary purpose of this thesis was to provide a 
method of checking the geometry and element connectivity 
input data for two finite element programs, ADINA and 
SAP IV. This preprocessor will accept the ADINA or SAP IV 
data deck, with minor modifications, and generate a graphi- 
cal display of the finite element model. The display is 
an oblique orthographic projection, and any orientation 
may be specified. Several options are available: exploded 
plots, partial plots, node numbering, element numbering, 
and others. Elements with three nodes on the same edge 
are plotted with a continuous curve on each edge generated 
by an interpolated parabola. Displacement postprocessing 


capability also exists. 
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I. INTRODUCTION 


A. GENERAL 


j 
The continued developments and advancements of the 
finite element method this last decade have provided 
greater machine capabilities than ever before. Here at 
the Naval Postgraduate School (NPS), two of the more ] 
| 
favorable finite element programs are SAP IV [reference 1] | 
for linear analysis and ADINA [reference 2] for nonlinear 
analysis. However, with the large amounts of numerical | 
input/output data and automatic mesh generation, it is | 
impractical to check and reduce this data without a graphi- 
cal representation. Data checking is divided into two 
areas: preprocessing and postprocessing. ‘ 
l. Preprocessing 
Preprocessing is the checking of the input data 
deck. Errors in a finite element program occur basically | 
in two areas. First, how close is the mathematical model | 
(boundary conditions, loading conditions, material proper- : 
ties, etc.) to the real problem? Second, are numerical 
errors present, or did misinterpretation of instructions 
occur in data deck preparation? Of the second type, most 
common errors are found in the geometry and element connec- 
tivity data. Preprocessing includes the forming of a 
graphical representation of the finite element model on 


which geometry and element errors are easily detected. When 
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node and element numbering options are available, the graph 
aids in the physical interpretation of the output data. 
Preprocessing is not a foolproof method of eliminating 
errors, but it does provide a tremendous advantage to the 
user. Preprocessors may be incorporated into the data check 
mode of the finite element program. However, to modify a 
large and complex program is dangerous. This may not be 
the best approach. A safer method is to develop a pre- 
processor which will read the finite element program deck 
separately with minimum modifications to that deck. 
2. Postprocessing 

Though not as important as preprocessing, post- 
processing is extremely helpful in output analysis. Probably 
the most common and useful type of postprocessing is the 
contour plot. Appendix C of reference 4 lists a program 
developed to produce contour plots of stress data from finite 
element models. Contour plots can easily be adapted to a 
2D system, but 3D requires plotting the contours on 2D 
surfaces, a bit more complicated. Two methods of post- 
processing of displacements are the plotting of a deformed 
model or placing scaled vectors at the nodes. In the cases 
where the displacements are small, multiplication by a 
magnification factor produces an exaggerated representation. 
Like preprocessing, postprocessing can be incorporated into 
the finite element program directly, or done separately. 
When done separately, the finite element program must still 


be modified slightly to obtain a punched deck of the stresses 


ll 


al 


and displacements in the desired format. This thesis is 


primarily concerned with preprocessing. 


B. HISTORY OF DEVELOPMENT OF PSAP1 

A package [reference 4] containing digital computer 
programs for generating oblique orthographic projections 
and contour plots was produced by the National Aeronautics 
and Space Administration's (NASA) Langley Research Center 
(LRC) and distributed by the National Technical Information 
Service in January, 1975. The programs are completely 
general. Both programs contain options for selecting 
various plotting equipment including CALCOMP, VARIAN, and 
cathode ray tube (CRT) displays. With minor modifications, 
they can be adapted to most any system. 

1. SUBROUTINE PSAP Implementation 

Losh [reference 6], for his master's thesis in 

aeronautical engineering, implemented the preprocessor and 
postprocessor program, PSAP, at NPS in December, 1976. 
Modifying the LRC package [reference 4], Losh adapted 
SUBROUTINE PSAP to the NPS IBM 360/67 system using the 
CALCOMP model 765 plotter. PSAP serves as a preprocessor 
for SAP IV models, and serves as a postprocessor for dis- 
placements of those models. Unfortunately, PSAP is severely 
limited in the type of elements it can plot. 

2. Motivation for SUBROUTINE PSAP1 


With the introduction of the ADINA [reference 2] 


program at NPS in January, 1977, and with expectation of 


doing future analysis on ceramic turbine blades, it was 


desired to expand PSAP to include all ADINA elements and 


the 8-20 node brick elements in SAP IV. Like PSAP, PSAP1 
contains preprocessing and displacement postprocessing 
capabilities. PSAPl is presented in this thesis and has 
the following improvements over PSAP: 

a. Preprocessing for all ADINA elements. 

b. SAP IV 8 and 8-20 node elements. 

c. Expansion of SUBROUTINE ERROR. 

d. Interpolation of curves using shape functions 
{reference 3] through three points on the edges of the 
8-20 node brick elements and the 4-8 node plane elements. 

e. Improvements in defining the plot origin. 

£. Addition of an option (ISCALE = 0) to plot 
sections of a model without losing perspective. 

g. Several other minor modifications considered 


improvements. 


C. PRESENT CAPABILITY 

PSAP1 has the capability to plot all ADINA elements 
and all SAP IV elements except the pipe element. It will 
interpolate curves on the edges of 4-8 node plane elements 
and 8-20 node brick elements. Many options are given in 
Appendix A. Some of the more frequently used options are 
listed below: 

1. Numbering of grid points (NOTAT = 1). 


2. Numbering of the elements (NOTAT = 2). 


3. Exploded plot (KDISP = 2). 


4. Postprocessing of displacements (NUDISP or NVDISP 


or NWDISP = 1) in two forms: plot of deformed structure 


(KDISP = 1) or displacements represented by vectors at 
the nodes (KDISP = 3; see reference 6). 

5. Symmetric representation about the XY (KSYMXY = l), 
XZ (KSYMXZ = 1) or YZ(KSYMYZ = 1) planes. 

6. Option to plot sections of the model (partial plot) 
to obtain a better view. Partial plots may be plotted to 
the scale of the complete model to avoid losing perspective 
(ISCALE = 0) or blown up to obtain a better view (ISCALE = l). 
Multiple plots may be obtained using the same geometry and 
same displacement data (KODE = 1), same geometry and new 
displacement data (KODE = 2), or new geometry and new 
displacement data (KODE = 3; see figure l). 

In general, multiple plots (sections, partial plots, 
additional problems) present no problem. Plotting package 
user courtesy dictates that no more than 5 plots be placed 
on the CALCOMP plotter at any one time. Also, if the plots 
contain many elements (especially 8-20 node elements), it 
is possible to run out of space in the plotting data sets. 


When this happens, you will receive 
ERROR IHC240I STAE, ABEND CODE IS: SYSTEM OB37 SYSPLOT. 


The best thing to do is split the run into two jobs. If 
the job must be run on one job (i.e., a large number of 
elements in the model or an assembly drawing where the 


scale of multiple plots is the same), then SYSPLOT space 


may be increased [references 7 and 8] by adding the card 
//GO.SYSPLOT DD SPACE=(CYL, (needed space) ) ,SYSOUT=C 

just prior to card 

//GO.FTLOFOO1l DD UNIT = SYSDA 


in Appendix A. It would be wise to seek advice from a 
consultant in Ingersoll 146 if additional plotting space 


is required. 


D. EASE OF MODIFICATION 

Both PSAP and PSAP1 are written to maintain as much 
generality as possible for ease of expansion and modifica- 
tion. Several FORTRAN statements, variables and subroutines 
are not used. They were left purposely unchanged. Although 
PSAP1 specifically reads ADINA and SAP IV data, it can 
easily be expanded to include any geometry and element data 
format. Simply study the read-in and storage methods (see 
Section II), and construct appropriate subroutines to read 


any particular format. 


II. PROGRAM ORGANIZATION AND DESCRIPTION OF OPERATION 


A. PSAP1 FLOW CHART 

Figure 1 is a condensed flow chart of PSAP1. Probably 
the most important information given on this chart is the 
sequence in which the data cards, NAMELIST OPTION and 
NAMELIST PICT are read. Remember, when generating a 
sequence of plots, once a parameter has been defined, it 
retains that value until it is reassigned. Note that when 
KODE = 1 or 2, the original values of NAMELIST OPTION and 
NAMELIST PICT are retained until they are changed. However, 
when KODE = 3 a new title card, NAMELIST OPTION and a set 
of problem data are read. All variables in NAMELIST OPTION 
and NAMELIST PICT are assigned their default values. A 
new problem begins in this case. It is important that the 
last NAMELIST PICT to be read must contain the value of 


KODE = 0. 


B. NAMELISTS AND EULER ANGLES 
1. NAMELIST OPTION 

Description and default values of NAMELIST OPTION 
are given in Appendix A. Basically NAMELIST OPTION varia- 
bles pertain to the given problem: the number of nodes, 
geometry format, displacement format, space between plots, 
and paper size. Default values for NAMELIST OPTION are 
set, and NAMELIST OPTION is read at the beginning of the 


problem. Assigned values will remain until exit from the 


Figure 1. Flow chart 
of program PSAP1 


Open plotting data sets. 
Define the origin. 


Read and plot title card. 

Assign defaults to &OPTION 
and &PICT. 

Read &OPTION. 

Read and print geometry 

and element data. 


KUDISP=KVDISP= 


Read and print displace- 
ment data. 


Read &PICT. 
Scale the plot. 

Rotate original coordinates. 
Plot the elements. 


Plot title at top of graph. 


Position the pen. 
Call PSTOP and close data sets. 
Stop 
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a 


program (KODE = 0) or a new set of problem data is read 
(KODE = 3). 
2. NAMELIST PICT 

Like NAMELIST OPTION, NAMELIST PICT variable des- 
criptions and default values are given in neeeadi« A. 
Basically NAMELIST PICT variables pertain to a given plot. 
One NAMELIST OPTION may apply to several successive plots, 
but each NAMELIST PICT defines a unique plot. That plot 
may include the whole model, part of the model and any 
options defined in NAMELIST PICT. NAMELIST PICT also speci- 
fies the viewing plane through the Euler angles (figure 18, 
Appendix A). 

a. Oblique Orthographic Projections (Euler Angles) 


An example of an oblique orthographic projection 


of a finite element model is given in figure 18 in Appendix 
A. The model can be viewed in any selected orientation. 
Euler angle transformations are used to specify orientation 
of the model to be projected. As described in reference 4, 
this transformation resolves the coordinate system of the 
model to a principal viewing plane (i.e., X_Y_, X_2Z2_, YZ.) 
on which the display is to be plotted. Prior to rotation, 
the model coordinate system (X,Y,Z) is coincident with the 
coordinate system containing the viewing planes (X57¥57%,)- 
The viewing planes are fixed, and the model is rotated about 
its model coordinate system. The rotations (,6,¢) of the 
body about the model axes (X,Y,Z) are shown in figure 18, 


Appendix A. The NAMELIST PICT variables KHORZ (horizontal 
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axis), KVERT (vertical axis), PSI (W), THETA (6), and 
PHI (9) specify the viewing plane and Euler angles. The 
order of the Euler angle rotations must be PSI, THETA 


and then PHI. Mathematical transformations are: 


Xo X 
Yo, = Gao fi tas 4 
26 Z 


[A] = Siny Cos w 0 
0 0 1 

Cos 8 0 Sin 6 

-Sin 8 0 Cos 8 
i 0 0 

[Ay] = 0 Cos -Sin $ 


0 Sino Cos 9 


SUBROUTINE ROTAT calculates the transformation matrices for 


every NAMELIST PICT, except if ISCALE = 0. ISCALE = 0 


directs the scale of the plot to be the same as that of the 


 salmemiann—-~--~-~-~-->~ eae een nna 


previous plot. Should a rotation occur with ISCALE = 0, 
the plot width could exceed the paper width. 
b. Scaling 

The safest scaling method is automatic scaling 
(ISCALE = 1). The user may specify a scale (ISCALE = 2) 
and the plot origin (XORGN,YORGN), but one must be careful 
not to run the plotting pen off the graph paper. ISCALE = 0 
is a very useful option. The plot will use the same scale 
as the previous plot. It is useful in an assembly graph 
where examination of a mesh in sections without losing 
perspective is desired. Example 3, Section III, illustrates 
the option ISCALE = 0. When ISCALE = 1 in a NAMELIST PICT 
defining a partial plot, a "blow-up" of that section is 
obtained. ISCALE cannot be zero in the first NAMELIST PICT. 

c. Partial Plots 

To develop a partial plot, three methods of 
segregating elements exist: first, by the X, Y, and Z 
cutting planes; second, by node numbers, and, third, by 
element numbers. If a model has an area where the elements 
are relatively small, a "blow-up" may be desired. Choose a 
numbering scheme or coordinates to define the section to 
be segregated using one of the methods above. Example 3 
(figure 14, Section III) uses X, Y, and Z cutting planes to 


define the partial plots. Example 4 (figure 17, Section III) 


uses element numbers to section the plots. 


C. NODAL POINT (GEOMETRY) INFORMATION READ-IN 


Nodal point data is read in by the GEOMn subroutines 
(GEOM1, GEOM2, and GEOM9, see figure 2). Since SAP IV 
and ADINA data decks are similar, SUBROUTINES GEOM1 and 
GEOM9 are also very similar. They are both constructed to 
read and generate data in exactly the same way as ADINA and 
SAP IV. Ail data not needed by PSAP1 is disregarded and 
the nodal point data is stored in array 222 (figure 3). 
After studying storage array ZZZ2 and GEOM1 (or GEOM9), 
a user moderately familiar with FORTRAN programming could 
easily construct a user supplied subroutine (GEOM2) to 


read the nodal point data in any desired format. 


D. ELEMENT (CONNECTIVITY) INFORMATION READ-IN 


After reading and storing the nodal point data, the 
element data must be read. The GEOMn subroutine will read 
the element control card (NPAR, references 1 and 2). 
SUBROUTINE ELTYPE (figure 2) calls the proper element sub- | 
routine to read the element data specified on the element 
control card. If several groups of elements are to be 


read, the process is repeated until all of the element 


groups have been read. Although the nodal point data is 
stored in array 22Z, the element connectivity is read and 


stored on device 10 (disk). 


E. DISPLACEMENT DATA READ-IN FOR DISPLACEMENT POSTPROCESSING 
Displacement data may be read in (figure 1) by SUBROUTINE 


DATA9 (KDATA=9) or SUBROUTINES DATA1 or DATAS (user supplied, 
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KDATA=1 or 5). When read, displacement data is stored 

in the last three columns of array 22Z (figure 3). PSAP1 

can postprocess displacements for both ADINA and SAP IV. 

The difficulty comes in obtaining a punched deck of cards. 
Reference 6 gives a description of how to obtain a deck of 
cards for SAP IV in a format acceptable to SUBROUTINE DATA9. 
ADINA has no such provision. However, when preprocessing, 

the displacement data will be omitted (NUDISP=NVDISP=NWDISP=0), 


and this step will be by-passed. 


F. PLOTTING LOGIC 

SUBROUTINE PLOTX (figure 4) is the main plotting routine. 
Since the nodal point data is stored in array 2Z2Z and the 
connectivity is stored on device 10, it is a simple matter 
to read the connectivity from device 10 (one element at a 
time), and connect the nodes as they are defined in references 
l and 2. For example, NEND = number of nodes defining the 
connectivity of a single element, NUMEL = the element number 
and NODE(NEND) is the array containing the connectivity. 
Device 10 contains this information successively for each 


element. The statement 


READ (10) NEND,NUMEL, (NODE (I) , I=1,NEND) 


will read the element connectivity to be plotted. The 
8-20 node brick (ADINA and SAP IV) and the 4-8 node plane 
elements (ADINA) may have 3 points defining each edge. If 


the midpoint node is defined, then isoparameteric 
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START 


Read elements individually to be 
plotted from device 10. 
READ(10) NEND,NUMEL, (NODE(I) ,I=1,NEND) 


Check if element is in bounds ” 
to be plotted 


YES 


NEND=2 NEND=4 NEND=8 NEND=12 OR 20 


Plot SAP IV 
8 node 
brick 
elements 


Plot ADINA 
8 node plane 
and 20 node 
brick elements 
and SAP IV 

20 node brick 
elements 


elements 


NO Check if all elements have been 
read from device 10. 


YES 


RETURN 


Figure 4. Flow chart for PSAPl plotting subroutine, PLOTX. 
Array NODE contains the connectivity of the 
element being plotted. NEND = the number of 
nodes defining the connectivity of a single 
element. NUMEL = the element number. 


shape functions (reference 3) are used to interpolate along 
the three-node edges. These shape functions are identical 

to those used by ADINA and SAP IV, so the geometry represented 
graphically is identical to the problem solved in ADINA and 
SAP IV. SUBROUTINE CURVE does the interpolation with the 


following equations: 


X =N1* X + N2 * X + N3 * X 


fo) 01 02 03 


Y =N1l*Y + N2 * Y +N3 * Y 


fe) ol 02 03 
Ni w.S * (S —- 1.0) 7 2.0 

WZ == (S + 1.0) * (8 = I20) 

NS = S * (S + 1.0) 7 2.0 


PSAP1 uses the NPS plotting package [reference 5] 


subroutines. 
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III. PSAP1 SAMPLE PROBLEMS 


The following examples have been chosen to illustrate 
some of the most useful options of PSAPl. Prior to 
attempting to use PSAP1, the user should have the problem 
defined and the cards punched in the format of references 
l and 2. Appendix A of this thesis gives a complete des- 
cription of deck preparation for PSAPl here at NPS. This 
section should prove helpful in the understanding and 


interpretation of the options presented in Appendix A. 


A. ADINA EXAMPLES 

1. Reinforced Concrete Beam (example 1, figure 5) 

This example was chosen because it illustrates 

the importance of the exploded plot when more than one 
element group is used. It is taken from the ADINA manual 
{reference 1]. Figure 6 is a listing of the data cards 
as they are prepared for ADINA. Figure 7 indicates how 
that deck would be modified for PSAP1. Note: load cards 
are removed, NAMELIST OPTION and NAMELIST PICT are added, 
and the title to be plotted on the graph is added in figure 
7. Otherwise, figures 6 and 7 are the same. Figure 8, 
part (a), illustrates an undistorted (KDISP = 0) PSAP1 plot 
with the nodes numbered (NOTAT = 1). Figure 8, part (b), 
shows the same mesh in an exploded form (KDISP = 2) with 
the elements numbered (NOTAT = 2). Note how the truss 


elements are visible in the exploded plot. 
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2. Flat Plate With Hole (example 2, figure 9) 


This is a well known problem with which one can 
calculate the stress concentration on a hole in a plate 


under axial tension. Figure 10 is a listing of the PSAP1 


deck set-up. The mesh is composed of ADINA variable 4-8 
node plane elements. Notice on figure ll, parts (a) and 
(b), how the interpolating shape functions round off the | 
3-node edges. Part (a) has the nodes numbered (NOTAT = 1). 
Part (b) has the elements numbereid (NOTAT = 2), and illus- 

trates the use of the symmetry option (KSYMXY=KSYMXZ=1). 

The symmetry option enables one to obtain a picture of the 


complete plate even though the model only consisted of a 


| 
quarter plate with proper boundary conditions. 
B. SAP IV EXAMPLES 
1. SAP IV Truss Problem (example 3, figure 12) 
| 
Figure 13 is a listing of the PSAP1 data deck. 
Figure 14 indicates how multiple partial plots can be used 
to obtain a better representation of the model. Part (a) 
of figure 14 is the complete model. Part (b) is the left 
half (XXMAX = 50'), and part (c) is the right half 
(XXMIN = 50', XXMAX = 1.0E20'). Notice also that for 
Parts (b) and (c), ISCALE = 0, which means succeeding plots 
have the same scale as the first. Had ISCALE equaled 1 
(blow-up), then the width of the half view would have been 


the same as the complete model. Figure 14 size is limited 


by the NAMELIST PICT variable PLOTSZ. 


Figure 9. Example 2, flat plate with a hole in 


tension. ADINA 4-8 


Thickness = 


Young's modulus 
Poisson's ratio = 
Total elements = 


Total nodes = 


node plane element. 
1 inch 

30.0 x 10° psi 
sas 

18 
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Distributed Load = 


1000.0 psi 
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Figure 11. Example 2, PSAP1 output graphs, page l of 2. 


42 43 as 


(a) Actual mesh, nodes numbered. 


‘ | 


Continued, page 2 of 2. 


Figure ll. 
(b) Symmetric representation, elements numbered. 
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Figure 14. Example 3, PSAP1 output graphs. 


(c) Right hand 
side. 


(b) Left hand 
side. 


(a) Complete structure. 


_ 2. Cylindrical Bar With Spherical Hole (example 4, 
figure 15) 


This mesh could have several uses. Two might be 
to calculate stress concentrations if the bar is under 
axial load or to calculate loading if the void is under 
pressure (i.e., dispersed nuclear fuel pellet). Figure 16 
is a listing of the PSAP1l data deck. Figure 17, part (a), 
is a representation of the complete model. Parts (b), 
(c), (d), and (e) of figure 17 are partial plots of the 
total structure using options of the undeformed structure 
(KDISP = 0) with node numbering (NOTAT = 1), and the 


exploded plot (KDISP = 2) with element numbering (NOTAT = 
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2). 


3"R Axisymmetric 


_ eylindrical bar. 


ae 


Z 
SA Spherical void under 
BOUNDARY CONDITIONS 12000 psi pressure. 
, = 
X fixed on YZ plane (xX 0) Young's modulus 6 


Y fixed on XZ plane (Y 30.0 X 10 psi 


Z fixed on XY plane (2 


“a wt 
oO 


Poisson's ratio = 
PERE 
Figure 15. Example 4, six inch diameter cylindrical bar with 
a two inch diameter spherical void on the center 
line under pressure. 
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Continued, page 6 of 6. 


Figure 16. 


- 
—_ 
YUN 
rr tm 
NWO 8 & 
wunoo 
ReeRavel t 
VaANMVEWO 
res ALO 
2002 4OW 
WSSYORKU 


Not part of input deck. 


kkEKK 


ssp aaa 


G 
Qu 
0 
4 
ion) 
v 
3 
Qy 
BS) 
3 
Oo 
qd 
Ay 
eg 
¥ 129) 
Ay 
SSE a 
~ aged 
wv 
i) 
qa 
Q, 
= ‘ 
© 
*% 
ica) 


at 


- ~ ~ ~ - - 


Figure 17. Continued, page 2 of 5. 


(b) Elements 1-10, nodes numbered. 


Figure 17. Continued, page 3 of 5. 


(c) Elements 1-10, elements numbered, exploded plot. 


Figure 17. Continued, page 4 of 5. 


(d) Elements 11-20, nodes numbered. 


Figure 17. Continued, page 5 of 5. 


(e) Elements 11-20, elements numbered, exploded plot. 


IV. CONCLUSIONS AND RECOMMENDATIONS 


With a little study of PSAPl and figure 2, one realizes 


how simple it would be to expand PSAP1 to include virtually 


an unlimited number of geometry data input formats. Simply 
supply the GEOMn subroutines, element subroutines and 

modify SUBROUTINE ELTYPE if necessary. Another interesting 
project would be to incorporate PSAP1 into ADINA and SAP IV 

so that a mesh plot could be obtained in the data check 

mode only. Reference 4 contains routines for plotting stress 
contours which could be incorporated. Hidden line logic 

is available. It would be difficult but could be incorporated. 
Hidden line plots are pretty, but they are not very practical 
when preprocessing because in preprocessing it is advantageous 


to see all of the nodes. 


APPENDIX A 


PSAP1 USERS MANUAL 


I. NAMELISTS 


It may be useful for the user to review NAMELISTs in any 
good FORTRAN manual. However, a short description of the 
NAMELIST input is given here. Only columns 2 through 80 
of each card will be read. The computer expects to find a 
special delimiter symbol in column 2 of the first card 
followed by the NAMELIST name. The delimiter symbol for 
the IBM 360-370 series is the ampersand & ; other processors 
use the dollar sign $. Following the first card comes the 
list of variables and their assigned values. Some important 
points about the variables are: 

A. Variables and their values need not appear in any 
special order. 

B. Predefined variables need not appear in the list 
if it is not necessary to change the value. 

C. A comma should follow each assigned value. 

D. It is recommended that each variable appear on a 
separate card to facilitate change. 

The last card following the NAMELIST variables contains, 
beginning in column 2, the delimiter sign followed by the 
word END. For example, &OPTION (read NAMELIST OPTION) 


contains 12 variables. All are initially assigned default 
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values. But, suppose three of those values (NNDEST=200, 


KGEOM=9, YSPACE=2.0) do not apply to a specific problem. 


The &OPTION would be read as follows: 


&OPTION 
~NNDEST=300, 
_KGEOM=1, 
_YSPACE=5.0, 
~&END 


Blank space in column l. 


The following sequence of cards is necessary to use 


PSAP1. 


NOTES 


(1) 


(2) 


(3) 


II. INPUT PROCEDURES 


IBM CARDS 
COLUMNS 


123456789-------------- 


// (Standard green job card] 
// EXEC FORTCLGP, REGION. GO=150K 


//FORT.SYSIN DD * 


Cc MAIN PROGRAM 
DIMENSION 222(NZ) ,DISPD(5,3,NON) 
CALL PSAP1(2Z2,NZ,DISPD,NON) 
STOP 


END 


/* 
//GINK.USDD DD UNIT=3330,VOL=SER=DISKO2, 
// DISP=SHR,DSN=S1153.PSAP1 
//UINK.SYSIN DD * 

INCLUDE USDD (LOADM) 


ENTRY MAIN 
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(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


/* 
//GO.FTLOFO01 DD UNIT=SYSDA, 
// SPACE=(CYL, (3,1)), 


// DCB= (RECFM=VS , BULKSIZE=3520) 


//GO.SYSIN DD * 


[PSAP1 title card as it appears on plot] 


&OPTION 


(&OPTION variables to be initially set or changed] 


& END 


{[ADINA or SAP IV geometry data. Title 
card through element data — remove load cards] 


[Case identification card — Omit if IDCASE = 0] 


[Displacement data cards — Omit for preprocessing 
only] 


&PICT 
{&PICT variables to be initially set or changed] 


&END 


{Additional data as defined by last value of 
KODE in NAMELIST PICT — Omit if last 
value of KODE = 0] 


/* 
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a tn ey EA rt Se et et ie tain: 


Notes: 


(1) Standard basic deck set up as described in chapter 
3 of reference 8. 


// (Standard green job card] 
// EXEC FORTCLGP, REGION.GO=150K 
//FORT.SYSIN DD * 


(2) Main program. 
DIMENSION Z22Z2Z(NZ) ,DISPD(5,3,NON) 
CALL PSAP1(ZZZ,NZ,DISPD,NON) 
STOP 
END 
The main program has two functions: to allocate 
fast storage space, and to call PSAP1l. NON must 
be greater than the number of nodes. NZ must be 
greater than 4*NON(7*NON if displacement data cards 


are to be input. i.e., NUDISP, NVDISP, or NWDISP = 1). 


(3) Using load module library. 
/* 
//LINK.USDD DD UNIT=3330,VOL=SER=DISK02, 
// DISP=SHR,DSN=S1153.PSAP1 
//GINK.SYSIN DD * 


INCLUDE USDD (LOADM) 
ENTRY MAIN 


PSAP1 should be precompiled and stored in the machine 


in a load module. Load modules are a type of user 


library described in reference 7 and chapter 3 of 
reference 8. Since PSAPl requires over one minute 


to compile, precompiling results in a large time 


| 
| 


64 


(4) 


(5) 


(6) 


Savings. These cards are subject to change, and the 
most current version is determined by the user main-~ 
taining this library. If PSAP1l is not on a load 
module, these control cards may be replaced by 
subroutine PSAP1 (which consists of a box and a 


half of cards). 


Allocation of storage. 
/* 
//GO.FTLOFOO1 DD UNIT=SYSDA, 
// SPACE=(CYL, (3,1)), 
// DCB=(RECFM=VS , BULKSIZE=3520) 
PSAP1 uses a slow storage device to store the element 


connectivity. These cards allocate 3 cylinders as 


described in chapter 3 of reference 8. 


Deck set up card. 


//GO.SYSIN DD * 


This is a standard deck set up card described in 


chapter 3 of reference 8. 


PSAP] title card. 


PSAP1 title card is the title that will appear on the 
graph. Make sure a user identification is on this 
card. It consists of 80 alphanumeric characters. 

The first 40 characters will form the first title 


line. The last 40 will form the second line. 
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(7) NAMELIST OPTION. 


(Note: start in second column) 


&OPTION 
[&OPTION variables to be initially set or changed] 
& END 
VARIABLE-DEFAULT DESCRIPTION 
NNDEST-200 Must be equal to the number of nodes. 
NUDISP-0 0 — X direction displacements not input. 
1 — X direction displacements input. 
NVDISP-0 0 — Y direction displacements not input. 
1 — Y direction displacements input. 
NWDISP-0 0 — Z direction displacements not input. 


1 — Z direction displacements input. 
(Note: unless displacement data is to be input, 
allow NUDISP, NVDISP, and NWDISP to default.) 
KGEOM-9 Specifies the geometry input format. 
1l- Subroutine GEOM1 reads in ADINA data 
deck geometry and connectivity. 
2- Subroutine GEOM2 may be supplied by 
user along with subroutines NSPLAN, NS3DEE and NSTRUS 
to read nodal data and connectivity in any format. 
9- Subroutine GEOM9 reads in SAP IV 
data deck geometry and connectivity. 
KDATA-9 Specifies the subroutine and corresponding 


method of input for displacement data. 


l- Subroutine DATA1, a user supplied 


subroutine. 


2- Subroutine DATA2, a user supplied 
subroutine. 
9- Subroutine DATA9, reads a punched 
output displacement deck from execution of SAP IV as 
presented in reference 6. (Note: unless displacement 


data is to be plotted, allow to default.) 


| NVALUS -0 Not incorporated, allow to default. 

| IRESEQ-1 Not incorporated, allow to default. 
KPLOT-1 Not incorporated, allow to default. 
YSPACE-2.0 Space between plots in the Y direction 


in inches when successive plots are plotted (i.e., 
KODE # 0). The graph title is plotted both on the 

top and at the bottom of each set of graphs controlled 
by a given NAMELIST OPTION. The space between the 
title and the plot is YSPACE/2.0. 

PSIZE-9.0 Paper size in the X direction, in 
inches. Used in scaling of the plots to insure this 
dimension is not exceeded. However, when manual 
scaling (ISCALE = 2; see NAMELIST PICT) this protection 
is not available, and it is possible to exceed the 
paper width. 

IDCASE-0 O- No identification card precedes 
the deck of displacement values. 

l- Identification card precedes the 


deck of displacement values. 


(8) 


(9) 


SAP IV or ADINA data cards. 


Here insert the geometry (node coordinates) and 
element connectivity. This includes the title card 
through the element data cards; the load cards are 
removed. For a SAP IV data deck, NAMELIST OPTION 
variable KGEOM = 9. For an ADINA data deck, 


KGEOM=1. 


From ADINA deck remove: 
1. Applied loads data. 


2. Frequency and mode shape calculations data. 


From SAP IV deck remove: 
1. Concentrated load mass data. 
2. Element load multipliers. 


3. Dynamic analysis cards. 


Otherwise, these cards are exactly the same as the 
deck prepared for SAP IV or ADINA. PSAPI1 is not 

limited to SAP IV and ADINA. The user may specify 
any unique format (i.e., KGEOM = 2). In this case, 


SUBROUTINE GEOM2 must be supplied by the user. 


Case ID card and displacement data cards. 


PSAP1 is intended to be used essentially for prepro- 


cessing. If it is desired to use the postprocessing 


option, the user is referred toreference 6. Otherwise, 


omit the case ID card and displacement data. 


(10) and (11) NAMELIST PICT. 


(Note: start in second column) 


&PICT 
[&PICT variables to be initially set or changed] 
&END 
VARIABLE-DEFAULT DESCRIPTION 
KHORZ-1 Integer designating the horizontal axes 
of the viewing plane. 1 = X,° 2= Yo: 3 = Zo: 
(See figure 18.) 
KVERT-2 Integer designating the vertical axes | 


of the viewing plane. 1 = Xo: 2 Se ae 

(See figure 18.) 
PHI~-0.0 Angular rotation of the model about its 

X axis in degrees (performed 3rd, see figure 18). 
THETA-0.0 Angular rotation of the model about its 

Y axis in degrees (performed 2nd, see figure 18). 
PSI~-0.0 Angular rotation of the model about its 

Z axis in degrees (performed lst, see figure 18). 
NEWFR-1 l- Frame change before plotting. 

2- No frame change before plotting. 

(Normally allow to default. A frame change resets 

the Y origin past the previous plot by YSPACE given 

in NAMELIST OPTION and resets the X origin at 0.0) 
ISCALE~1 O- No scale change. Use the same scale 


as the previous plot. Useful in an assembly graph 


where it is desired to examine a mesh in sections 


(°LZT ebed ‘pf soueTSJeA WOAF usyeL) 
-suetTd BburmetaA Z-x uT umoys uoTzOefoIrd OTydeszboyzI0 enbtTqGo 
ue ZOJ (SUOT}e}0A) SoTHue AeTNy pue sweqzsAs ajZeuTpI00D “gt seaNnbTYy 


% 


without losing perspective. See example 3 in 
Section III of this thesis. ISCALE cannot be zero 
in the first NAMELIST PICT. 

l~ Automatic scaling of plot and 
computation of proper origin location. 


2~ User specified origin and scaling. 


PLOTSZ-10.0 Maximum dimension desired on com- 
pleted plot prior to rotation. After rotation it 
is possible for the maximum dimension in the pro- 
jected plane to exceed PLOTSZ. The maximum size 
of the projection is unlimited on the vertical axis 
but is limited by the paper width (PSIZE in NAMELIST 
OPTION) on the horizontal axis. If rotation of the 
model causes the projection to go off the paper, 
it is rescaled prior to plotting. PLOTSZ is used 


in scaling only if ISCALE = 1). 


XORGN-0.0 X location of the plot origin. 


(Used only if ISCALE = 2.) 


YORGN-0.0 Y location of the plot origin. 


(Used only if ISCALE = 2.) 


PSCALE-1.0 Model size reduction factor. 
PSCALE is equal to the actual model size divided by 
the desired plot size. It is used only if ISCALE = 2. 
(Note that when ISCALE = 2 is used, no rescaling 


occurs if rotation causes the projection to exceed 


the paper width.) 


NOTAT-0 O- No numbering on plots. 


l- Numbering of grid points. 


2- Numbering of elements. 


XLHT-0.14 Height in inches of the integers 
specified by NOTAT. It should be a multiple of 
-07. If XLHT is not a multiple of .07, it will be 
rounded to the nearest multiple. XLHT has a 


maximum of .49 and a minimum of .07. 


KDISP-0 O- Plot of undeformed structure. 
l- Plot of deformed structure. 


2 


Exploded plot. 

3- Displacement represented by 
vectors. 

KDISP = 1 or 3 represents a form 
of postprocessing and displacement data must be 
input in (9). If postprocessing is desired, refer 


to reference 6. 


IDMAG-2 l- Direct magnification of displace- 
ment data by DMAGS. 
2- Scaling of displacement data 


to a maximum value of DMAGS. 


DMAGS-1.0 Magnification of displacements 
(if KDISP = 1 or 3). Reduction factor of elements 
(if KDISP = 2). 

KSYMXY-0 l- Symmetry about X-Y plane. 
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KSYMXZ-0 l- Symmetry about X-Z plane. 


KSYMYZ-0 l- Symmetry about Y-Z plane. 
A plate quadrant with KSYMXZ and 
KSYMXY equal to 1 would yield a complete plate. 


See example 2 in Section III in this thesis. 


To develop a partial plot, three methods of 
segregating elements exist: First, by X, Y and Z 
cutting planes; second, by node numbers; and third, 
by element numbers. The next ten variables are used 


to separate elements into partial plots. 


XXMAX , YYMAX, ZZMAX-1.0E20 Local cutting planes 


XXMIN, YYMIN, ZZMIN- (-1.0E20) Parallel to the principal 
planes. 


NDMAX-9999999 Maximum gridpoint identification 


number to be included in the plot. 


NDMIN-O Minimum gridpoint identification 


number to be included in the plot. 


NELMAX-9999999 Maximum element number to be 


included in the plot. 


NELMIN-0O Minimum element number to be 


included in the plot. 


| 
| 


KODE-0 


Specifies the control option after 
the plot is completed. 

O- Last plot, exit from program. 

l- Read another NAMELIST PICT. 

2- Read a new set of displacement 
data. (Postprocessing only.) For.KODE = 2, dis- 
placement data must be followed by another NAMELIST 
PICT. 

3- Read a complete new set of input 
data starting with a title card. 

For KODE = 1, 2 or 3, additional 
sections of the deck must be repeated. The deck 
must end with a NAMELIST PICT having a value of 


KODE = 0. 


Note: A most important point to remember when 
generating a sequence of plots is that once a 
parameter has been defined, it retains that value 
until it is reassigned. For example, if PLOTSZ is 
assigned a value of 8.0 for the first of a series 
of plots, and it is not redefined in any subsequent 
NAMELIST PICT; the value of PLOTSZ will be retained 
as originally specified. However, when KODE = 3 
and a new title, NAMELIST OPTION and problem data 
are read, all variables in NAMELIST OPTION and 


NAMELIST PICT are assigned their default values. The 
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aie RE A ers 


problem starts over in this case. Refer to the 


flow chart in figure 1 and study the path for 


different values of KODE. 


Delimiter card. 


/* 


Delimiter card is defined in chapter 3 of reference 


8. 


TII. SPECIAL FEATURES OF PSAP1 


A. POSTPROCESSING 

Reference 6 contains information and examples on the 
use of the postprocessor. SAP IV has the capability to 
punch displacement data cards in an acceptable format for 
PSAP1 (specifically subroutine DATA9). As of this writing, 
ADINA does not have this capability, and data cards would 


have to be punched manually by the user. 


B. PARTIAL DATA 

Reference 6, on page 30, establishes a procedure by 
which it is possible to input only a portion of the finite 
element model for a data check. This feature is valuable 
in a case where several different people are preparing 
different parts of a large data base for a problem and 
desire to check individually their inputs graphically for 
accuracy. All nodal coordinates for the entire model may 
be input, or only those that specifically define the portion 
of the finite element model to be plotted. In either case, 
all nodal coordinates that define the elements to be plotted 
must be specified. To use this feature the "element control 
cards" (described in references 1 and 2) must be modified. 
For example, if only the connectivity for elements 15 
through 50 of a problem are available, the changes below 


would be made. 
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1. All SAP IV Elements 


All SAP IV element control card changes would 
have similar changes because columns 6-10 contain the 
total number of group elements. Columns 66-70 are not 
used. Thus to plot only elements 15 through 50, make the 
following changes: 

a. Enter the upper bound (i.e., 50) in columns 
6-10. 

b. Enter the lower bound (i.e., 15) in columns 
66-70. 

2. ADINA Truss, 2D and 3D Elements 

ADINA element control cards for the truss, 2D 
continuum and the 3D continuum elements would be changed 
as follows: 


a. Enter the upper bound (i.e., 50) in columns 


b. Enter the lower bound (i.e., 15) in columns 
53-56. 
3. ADINA Beam Element 
ADINA element control cards for the beam element 
would be changed as follows: 


a. Enter the upper bound (i.e., 50) in columns 


b. Enter the lower bound (i.e., 15) in columns 


65-68. 
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